Field evolution of vortex lattice in LuNi2B2C seen by decoration in fields up to 1.5 kOe 
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Evolution of the vortex lattice, the rhombus-square transition included, for LuNi2B2C in the 
field along the c crystal axis, is tracked by the decoration technique pushed up to the record high 
(for this method) field of 1480 Oe. Decoration images are analyzed with the help of the Fourier 
transform, which indicates disordered structures in small fields of a few Oe. In fields H < 200 Oe the 
coexisting domains of different structures are observed. A technique based on the Fourier transform 
is employed to see the domains separately. The transition to the square lattice is recorded near 900 
Oe. The results are in agreement with predictions of the nonlocal London theory. 

PACS numbers; 74.60.-w, 74.60.Ec, 74.60. Ge 



I. INTRODUCTION 



The family of the rare-earth (i?) metallic nickel- 
borocarbides RNi2B2C has been the subject of consider- 
able recent attention due to the interplay of magnetic and 
superconducting properties. Bel The critical temperature 
Tc of these materials depends on the particular rare earth 
(or their combination), or on material purity, and ranges 
from zero to about 16 K in nonmagnetic members such as 
LuNi2B2C. The vortex lattices (VL) in these compounds 
exhibit a variety of phases with changing magnetic field 
and temperature. Studies of VL's are facilitated by avail- 
ability of large high quality single crystalsci with a rela- 
tively large Ginzburg-Landau parameter, k 10 — 20, 
low pinning, and a broad region of the HT phase dia- 
gram where the superconducting magnetic properties are 
nearly reversible. 

The structural phase transition from the triangular to 
square VL driven by the magnetic field is of a particular 
interest. A variety of experimental techniques are used 
to study tliiS|-|transition: small-angle neutron scatter- 
ing (SANS)ErEl, scanning tunneling microscopy (SXIV[)eI, 
and Bitter decoration by small magnetic particles .qLI The 
SANS and STM methods require rather high magnetic 
fields (few kOe and higher), whereas decoration is com- 
monly employed for fields on the order of several hundred 
Oe or less because of a relatively low resolution. 

For ErNi2B2C, the phase transition from the triangu- 
lar to square VL in increasing field has been recorded 
by SANS and complemented by decoration.0 Decoration 
patterns revealed coexistence of the square and triangu- 
lar lattices in fields near 600 Oe. In LuNi2B2C crystals, 
the square vortex lattice for hijeh magnetic fields of a 
few Tesla has been seen in STM.El No experiments in low 
magnetic fields were carried out for this material. The 
theoretical model, which describes most of the features 
of evolving VL's is based on the London theory corrected 
for nonlocal effects E 

In this work the VL's in LuNi2B2C were studied by the 



decoration method. To cover the field range of greatest 
interest for the VL evolution, the method was extended 
up to fields of 1.5 kOe, the field region unprecedented for 
the decoration technique in materials with k ^ 1. To our 
knowledge, the maximum field at which_the decoration 
has been attempted was 1100 Oe on Nb.EI 



II. EXPERIMENTAL 



We have used single crystals LuNi2B2C of the size 
~ 3x5 mm in the a6-plane and ^ 0.5 mm thick. 
The crystals-. were grown from Ni2B flux, as described 
elsewhere. cniJ Experiments were carried out in the field 
cooled regime with the field nominally parallel to the c 
axis (within a few degrees). For all fields except 1480 Oe, 
the temperature was set 4.2 K initially but may have 
risen to 7-9 K during the decoration experiment. For 
1480 Oe, the initial temperature was 1.5 K. To minimize 
effects of inhomogeneities and pinning, decorations were 
done on a particular location of the same sample which 
was cleared by acetone after each decoration. A scan- 
ning electron microscope was used to resolve decoration 
patterns on the sample surface and to determine VL ori- 
entations relative to the crystal by comparing decorations 
with electron channeling patterns (ECP). 

Routinely, the decoration experiments are done in 
fields near the lower critical field Hci on samples as thin 
platelets in the perpendicular magnetic field. Since there 
is no quantitative theory to estimate the resolution limit, 
we use the following qualitative guidelines in pushing the 
method to higher fields: 

1. The size of magnetic particles should be small rel- 
ative to the London penetration depth A(r) and to the 
intervortex spacing ^ y^(j)o/B. 

2. The density of particles should be small to avoid 
clustering. 

3. The field gradients due to VL's should be large 
enough for the interaction with the particle magnetic mo- 
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ment n to exceed the thermal energy, fj.H > 3k bT. 

The first condition is reahzed by manipulation of the 
buffer helium pressure within w O.ltorr. The optimal 
size of particles of 5 — 10 nm should not limit the resolu- 
tion up to fields 10 kOe. The second requirement is met 
by empirically minimizing the quantity of the evaporated 
magnetic material and by varying the distance between 
the heater and the sample. The most difficult to control 
and to satisfy is the last condition because the magnetic 
moment fi of iron particles is unknown. Besides, the ac- 
tual temperature of decoration is uncertain because of 
the heating during the iron evaporation. Still, one has to 
keep this T as low as possible. To achieve this, we used 
onset temperatures 1.5 — 2 K for decoration in high fields. 
In this work, we have found that the contrast images of 
decorated VL's in LuNi2B2C can be readily resolved in 
fields up to 1480 Oe. 

The information one can extract from decoration ex- 
periments in fields below Hd differs from the case H > 
Hci- In the first case, the intervortex distance exceeds 
the penetration depth A, and one can well resolve single 
vortices. However, the intervortex interaction in this sit- 
uation is too weak to arrange vortices in a well ordered 
lattice, given the competition with disordering pinning. 
Nevertheless, one can utilize the decoration data to ob- 
tain a rough estimate of the penetration depth. The field 
distribution of a single vortex in the bulk has a character- 
istic diameter oiJLX. Near the sample surface the vortex 
field "opens up" .Ej As a result, the field gradients at the 
sample surface are lower than in the bulk. The field lines 
at a distance A from the vortex core in the bulk, cross 
the surface at a larger distance from the core. One may 
consider the magnetic diameter of the vortex at the sur- 
face as L fa mX, where the empirical expansion factor m 
can be determined measuring the size L of the iron par- 
ticle clusters forming at each vortex site during the low 
field decoration. Data show that to « 4 provides a good 
estimate of A. 

In applied fields well above Hd, the fields of neigh- 
boring vortices overlap, and the resolution of decora- 
tion experiments drops. Nevertheless, the contrast of 
the decoration patterns still provides enough quantitative 
information on the VL symmetry and intervortex spac- 
ing. The computer processing of the images (fast Fourier 
transform and image filtering followed by inverse Fourier 
transform) greatly enhance our ability to study VL's and 
the domain structures in fields which were thought inac- 
cessible to the decoration technique. In this work, VL's 
were successfully decorated and analyzed in a broad field 
interval from 4 to 1480 Oe. 



III. RESULTS 

As is seen in Fig. |], in small fields vortices do not form 
a well ordered lattice. A straightforward count shows 
that the average intervortex spacing corresponds to one 



flux quantum per vortex. 

High magnification images of single vortices reveal the 
average size of the iron clusters of about 600 nm. Tak- 
ing this as an estimate of 4A at T w 0.5 Tc, we ob- 
tain Ao ~ 130 nm utilizing the two-fluid approximation 
A2(r) = Ag/(1 - T^/T*). Tius value is larger than 71 nm 
obtained from the Hd datajlj and ~ 100 nm (Ref. |l^) ob- 
tained from the magnetization measurements, but close 
to the estimate of Ref. |^. 




FIG. 1. Vortices in field H — 4.5 Oe along the c crystal 
axis. The Fourier transform in the inset indicates that the 
vortex system is an amorphous solid. 




FIG. 2. Vortex system for H = 40 Oe. Domains of well or- 
dered lattices are seen. The maxima of the Fourier transform 
are smeared due to varying orientations of the domains. 

In higher fields, intervortex interaction is stronger and 
vortices form ordered lattices. Starting from «40 Oe and 
up to «200 Oe we observe domains of nearly regular tri- 
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angular VL's; see Fig. VL orientations within each 
domain are quite random. Typical size of a domain is a 
few intervortex distances and increases with the magnetic 
field. 



An example of VL at H = 200 Oe is shown in the 
upper panel of Fig. ^. The Fourier transform (FT) of 
this pattern shows that the structure consists of three 
nearly regular triangular domains with different VL ori- 
entations. 




FIG. 3. The upper panel: the three-domain structure for 
H = 200 Oe. The domains (marked as 1,2,3) are shown sep- 
arately in Figs. ^. The lower panel: Fourier transform of 
the upper panel pattern indicates the presence of three well 
ordered domains. Each of the domains is responsible for 6 
maxima forming nearly a hexagon. 



By choosing the FT maxima which belong to one of 
the domains and clearing all other reflexes one can ob- 
tain with the help of the inverse FT the patterns of all 
three domains in real space separately. These patterns 
are shown in Fig. 0. 




FIG. 4. The three panels show three different VL's. Each 
one is obtained by the inverse Fourier transform of a selected 
sets of six maxima out of 18 of Fig. |^. 

Note that the domain size (in units of intervortex dis- 
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tance) for 200 Oe is large as compared to domains in the 
field of 40 Oe. 

The influence of crystal inhomogeneities and 
anisotropy on the VL domain structure., and ordering 
has been extensively studied for NbSe2.E3o The nonlo- 
cality brings a new complexity to the process of domain 
formation because VL^ now depend on T, H and the 
local mean-free path.Ej Comparing the ECP picture of 
Fig. ^ with Figs. U, we find that one of the closepacked 
directions on the upper panel is aligned with (110) direc- 
tion of the crystal, whereas that of the middle panel is 
aligned with (110). It is worth noting that the model of 
Ref. H predicts that the short diagonal of the rhombic VL 
cell in small fields should point in one of these directions. 




FIG. 5. Electron channeling pattern for the same place as 
the decoration of Fig. |^. 

In fields more then 200 Oe, the typical domain size con- 
siderably exceeds the intervortex distance. The orienta- 
tion of the closepacked VL direction in most domains cor- 
relates (within few degrees) with (100) crystallographic 
direction; this is seen in Fig. |^ in which the decoration 
pattern for H = 300 Oe and the ECP picture are pro- 
vided. 

Still, in other domains the closepacked direction co- 
incides with (110), the fraction of which decreases with 
field and aX H = 600 Oe reaches « 1/3. The VL unit cell 
is nearly rhombic with the apex angle (3 exceeding 60°, 
as can be extracted from the FT in Fig. The angle (3 
increases with field. The preferable orientation of the VL 
corresponds to the unit cell diagonal opposite to (3 being 
aligned with (100). 

In fields of about 600 Oe and higher, in some surface 
locations we observe domains of the square VL, whereas 
majority of the domains contain strongly distorted trian- 
gular (or, better to say, rhombic) VL, see FT in Fig. |^. 
The fraction of domains with the square lattice increases 
with increasing field. A similar situation has been seen 
in SANS data on VsSi, where the presence of domaiaa 
was deduced from analysis of the scattering maxima.ta 
Here, we see the domains directly. 




FIG. 6. The upper panel shows the decoration pattern for 
H — 300 Oe with the Fourier transform in the inset. The 
lower panel shows the corresponding electron channeling pat- 
tern indicating that one of the VL close-packed directions 
coincides with (100) crystal direction. 




FIG. 7. The applied field is H = 710 Oe. The direction 
analysis of the Fourier transform in the inset gives l3 ~ 82° 
and the corresponding diagonal of the rhombic cell aligned 
with (100). 
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The square domains are dominant in fields of about 
1000 Oe and higher, Fig. ||. The square diagonals are 
parallel to (100) and (010) crystallographic direction. 




FIG. 8. The vortex lattice in the LuNi2B2C single crystal 
for the field 1480 Oe parallel to the c axis. 

The results of our study of VL's in LuNi2B2C for a 
two particular locations on the crystal surface are col- 
lected in Fig. ^ where the field dependence of the angle 
j3 is shown. The triangles and circles show the decoration 
results from two different locations of the crystal. Few 
points for the same field indicate the scattering of the an- 
gle /3 within the decoration area of about 100 x 100 /im. 
In the field interval 200-600 Oe, the apex angle (3 of the 
rhombic unit cell slowly increases with field. With the 
further field increase, VL's deviate even faster from the 
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FIG. 9. The apex angle (3 versus H. The triangles and the 
circles are obtained at two different locations on the sample 
surface. A few triangles (circles) for the same field indicate 
that the data points are extracted from a few different loca- 
tions of the same decoration patch. 



IV. DISCUSSION 



The theory of Ref. || predicts that in fields less than 
Hi ^ 200 Oe along (001) crystal direction, the VL should 
have a rhombic unit cell with a short diagonal along ei- 
ther (110) or (110). The apex angle (3 of the rhombus 
should be less than 60°. Experimentally, up to fields 
about 200 Oe we did not record a well ordered VL. In- 
stead, a structure made of many domains with different 
lattice orientations is observed. The prediction for the 
value of /3 < 60° is also hard to verify in this field region. 

One can argue that in small fields energy differences 
between different possible VL configurations are compa- 
rable with characteristic pinning energies. In this situa- 
tion, the VL structure is affected by distribution of pin- 
ning sites which may vary throughout the sample. There- 
fore, we need better crystals of LuNi2B2C to support or 
disprove the low field theoretical scenario for equilibrium 
VL's for this material (for YNi2B2C this scenario has 
been confirmedO by SANS data in fields up to 1200 Oe). 

Still, the low field decoration results can be used to 
estimate the London penetration depth A at decoration 
temperatures of (0.3 — 0.5)rc. The value of A obtained by 
measuring sizes of the particle clusters, is in a reasonable 
agreement with values obtained by other techniques E'EZI 
It is to be noted that A estimated by measuring the "vor- 
tex diameter" d in decoration experimerLts,, is in a fair 
agreement with other data also for YBC0,t3t3 for which 
an empiric relation d = AX yields Aq ~ 150nnjj_and for 
NbSe2 for which the procedure gives ~ 270 nm.cd 

With the magnetic field increasing, the role of pinning 
diminishes because the intervortex interactions become 
stronger. The field region with H exceeding approxi- 
mately 300 Oe is more amenable for checking the theoret- 
ical model of Ref. ||. The model is based on the nonlocal 
London approach and as such requires a number of mi- 
croscopic input parameters (the Fermi surface, the pen- 
etration depth, the mean-free path, to name a few) for a 
quantitative comparison with experiment. Here, we use a 
simpler "mean-field" approach based on the assumption 
that the structural transition from the rhombic to the 
square VL at the transition field we call H2 is of the sec- 
ond order. This assumption is supported by numerical 
simulations of Ref. || which show that the rhombus-to- 
square transition at H2 is accomplished by an infinites- 
imal deformation so that at the transition the structure 
changes continuously. Moreover, the square phase has an 
extra symmetry element: the 4-fold rotation axis which 
is absent in the rhombic cell. Both of these features are 
characteristic of a 2nd order transition. We, therefore, 
can apply the general Landau theory according to which 
the energy density near the transition at H2 can be writ- 
ten as 



F{77; B} = Fo{B) + a{B - H2W + Wl'^ ■ 



(1) 
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Here, Fg is the part unrelated to the VL structure; the 
magnetic induction B for platelet-Uke samples is equal 
to the perpendicular applied field H. The field indepen- 
dent and positive coefficients a, b are related to a partic- 
ular VL; 77 is a properly defined order parameter which 
should be zero in equilibrium above H2. Since the apex 
angle of the high- field structure is tt/2, one can choose 
?7 = 7r/2 — /3. Note that the field B in (Q) plays the role 
of temperature in the standard Landau theory. 

The equilibrium structure corresponds to dF/drj =^ 
which yields = above H2 and 

for B < H2- The two signs correspond to two possible 
structures with apex angles /?^^^ -I- /S^^^ = tt; these are 
rhombic VL's rotated relative to each other by 7r/2. We 
then expect (7r/2 — /?) oc ^/H2 — B near i?2- In fact, the 
numerical simulations show that with good accuracy this 
dependence holds all the way down to small fields where 
P should approach 7r/3, i.e., ^/aH2/b w 7r/6. Thus, the 
apex angle (3 (in degrees) should satisfy an approximate 
relation 

/3 « 90° - 30 Vl - ^/^2 , B<H2, (3) 
/3 = 90° , B> H2 

in a broad field domain near H2. This formula is con- 
venient for comparison with data since it contains only 
one fitting parameter H2. The solid line in Fig. ^ is 
calculated with the help of Eq. (||) with H2 = 930 Oe. 

There are few reasons for a considerable spread of /3 
values shown in Fig. ^. Since the lattice should undergo 
restructuring during the cooling, the cooling rate may af- 
fect the resulting structure. The gas pressure may influ- 
ence the actual decoration temperature since during the 
"evaporation" of iron particles the sample heats up and 
cools down with the pressure dependent rate. Clearly, 
the spread related to different loci is related to crystal 
inhomoeeneities. Moreuver, the nonlocal London model 
predicta3 and the dataEJ confirm that the transition field 
H2 is sensitive to the crystal purity {via the mean- free 
path £). Even weak inhomogeneities in £ within one 
decoration patch of 100 x 100 firn^ may influence VL 
structures. As is seen from Eq. (|^), d[3/dB diverges as 
B — !■ H2, in other words, a small spread in values of H2 
may cause a considerable variation in [3. This statement 
can also be formulated as enhanced "structural fluctua- 
tions" near the second order transition at H2. In fact, our 
data show that the (3 spread peaks in the field interval 
700 - 1000 Oe, i.e., just near H2. 
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